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classification, the limits are geometrically motivated in terms of transformations of roots of 
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1 Introduction 

The discovery and classification of superintegrable systems is currently undergoing significant 
activity, see [22] and citations within. Superintegrable systems are dynamical systems, for 
example classical or quantum Hamiltonian systems, that have more conserved quantities than 
degrees of freedom. One of the many active areas of research in superintegrable systems is in their 
classification. On this front, there has been significant progress on maximally superintegrable 
systems in 2 and 3-D complex Riemannian manifolds with constants of the motion that are 
at most second-order in the momenta, a necessary condition for multi-separability. In two 
dimensions, all such systems have been classified [12]. Somewhat surprisingly, all such systems 
can be obtained through singular limits of a “generic” system defined on the sphere. This fact 
was first recognized by Bocher [1] in his search for metrics that admit the maximal number of 
second-order Killing tensors. Again in the 2D case, the limits between such systems were recently 
given explicitly in [19], where they were shown to induce contractions on the associated symmetry 
algebras as well as limits in the representations of such algebras. These representations are given 
in terms of classical hypergeometric orthogonal polynomials and the limits of the superintegrable 
systems correspond to the limits within the Askey scheme. 

The aim of this paper is to give the corresponding limits of superintegrable systems in 3D. 
We will show explicitly that all such systems are limiting cases of a “generic” system defined 
on § 3 and give the coordinate transformation that generates each limit. In order to describe 
the limits we discuss the classification theory of 3D second-order superintegrable systems [2, 3]. 
The classification makes use of classical invariant theory and associates to each superintegrable 
system a 7-dimensional representation of the conformal group. By studying the action of the 
conformal group on this representation, we can identify 10 equivalence classes each of which 
corresponds to a class of superintegrable systems. From the invariant classification, the necessary 
limits are abstractly motivated as opposed to the constructions of [19], which are given ad hoc. 
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The structure of the paper is as follows. We begin in Section 2 with an overview of the results 
in [2, 3] giving the invariant classification of nondegenerate second-order superintegrable systems 
in 3D. Section 3 contains the results of the paper; namely we prove that all such superintegrable 
systems are limits of a “generic” system by demonstrating the appropriate limits. Section 4 
contains a brief discussion of the results and future applications of this research. 

2 Invariant classification of superintegrable systems 

In this section, we describe the invariant classification of nondegenerate second-order superin¬ 
tegrable systems in 3D, where the metric is assumed to be conformally flat. The dynamics are 
given by the following classical Hamiltonian on 6-dinrensional phase space (p, x) 

3 

ft(p>x) = Y 9 lJ (^)piPj + V'(x). 

i,j =1 


The conformally flat assumption implies that the metric tensor takes the form g tJ = A(x)h*. 
Any function on the phase space A(p, x) will transform in time as 


dA 

dt 


{A,U} + 


dA 
~dt ’ 


where { , } is the standard Poisson bracket 


3 


{A,B} = Y 

i —1 


dA dB_ 

dxi dpi 


dB_dA 

dxi dpi' 


Thus, a function on phase-space will be invariant under the Hamiltonian flow if and only if 
it Poisson commutes with PL, ({A,PL} = 0) and is referred to as a constant of the motion. 
Conversely, the Hamiltonian is also invariant under flow generated by the constant of the motion 
and so the function A is also often referred to as a symmetry of the Hamiltonian. 

A Hamiltonian on 2n-dimensional phase space is said to be Liouville integrable if there ex¬ 
ist n mutually commuting constants of the motion that are functionally independent. We say 
that a Hamiltonian is superintegrable if there exist more than n constants of the motion and 
maximally superintegrable if there exist 2n — 1 such symmetries. Not all of the symmetries can 
mutually commute and so we drop the requirement that the integrals commute among them¬ 
selves. There is a close connection between superintegrability and noncommutative integrability, 
notably that maximally superintegrable systems are noncommutatively integrable. Indeed, the 
celebrated result of Nekhoroshev [23] concerning foliations of phase space for noncommutative 
integrable systems includes, as a special case, the result that, generically, bounded trajectories 
of superintegrable systems are closed. 

We would like to classify Hamiltonians PL on a conformally flat 3D manifold with 5 con¬ 
stants of motion that are second-order in the momentum, i.e., that there exist 5 (including the 
Hamiltonian) functionally independent constants of the motion 


3 

A = Y alj ( x )PiPj + IK(x), a lj = a?\ 
i,j=l 


such that 


{A,H} = 0. 
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If such symmetries exist then a direct computation of the Poisson commutators leads to the 
following set of differential equations: the Killing equations, which require the second-order 
part of A to be a Killing tensor on the manifold, and the Bertrand-Darboux (BD) equations. 
The BD equations are the compatibility conditions for the differential equations for W (x) and 
comprises a set of four linear second-order PDEs for the potential K(x), one for each integral 
that is not the Hamiltonian. Combining the equations into vector form gives 


B(a ij 


(V 22 - Kn\ 
D33 - 

Dl2 

Kl3 

V ^23 


( V A 

C(X,X ! k,a l 3 ,a l j e ) I V 2 \ , 


where the subscript denotes the partial derivative with respect to Xk and B and C are 
matrix functions of the given arguments of dimensions 12 x 5 and 12 x 3, respectively. Suppose 
in addition that the constants of the motion are functionally, linearly independent. This implies 
that the BD equations are of rank 5 and it is possible to solve for second-order derivatives of V 
as in 


( v ,u\ 


n 

-AS 1 - Rf - 

R 13 

2 S 2 + R\ 2 

2 S 3 + Rf 

\ 

v 22 


i 

2S l + Rf 


-AS 2 - R\ 2 - Rf 

2 S 3 + Rf 


V, 33 


i 

2S l + Rf 


2 S 2 + Rf 

-AS 3 - Rf - 

Rf 

v V2 


0 

Rf - 3S 2 


Rf - 35 1 

Q 123 


V.V3 


0 

R{ 3 - 3 S 3 


Q 123 

Rf - 3 S 1 


V^23/ 


\0 

Q 123 


Rf - 3 S 3 

Rf - 3S 2 

/ 


fV,ee\ 


Kee = ^(Kll 





V = 

V tl 

1 

+ v 22 

+ V 33 ). 




\vj 








Here V ee is not a second derivative, but a symmetry adapted variable. The use of V ee introduces 
redundancy in the equations (6 equations instead of 5), which is preserved for the sake of 
symmetry in the variables. If the compatibility conditions of (1) are assumed to hold identically 
then the potential depends on 5 parameters, the values of V ee , V\, V 2 , V 3 and V at a generic 
point, the last point coming from a trivial additive parameter. Such a system is said to be 
non-degenerate. For a non-degenerate system, the potential is uniquely determined by the value 
of the matrix in (1) at a given generic point xq. The values of K(xo), V] (xq), K 2 (xq), V^xo) 
and Vji (xq) are the parameters in the potential. As is clear from the form of the matrix, it 
depends on 10 functions 


f/") 123 q 1 q2 q 3 pl2 pl2 pl3 p 13 p23 p23\ — rn o T3 \ 

\ W , O , O , O , Ki , it 2 ,-ai , #3 , , -^3 / = {* 4 , -tt), 


which are assumed to satisfy the compatibility conditions for (1) identically. The compatibility 
conditions of (1) are first-order equations in the variables whose compatibility conditions are 
themselves identically satisfied. Therefore, the values of functions {Q, S, R} at a generic point xo 
are enough to uniquely determine the functions themselves along with a superintegrable system. 


Example 1. The harmonic oscillator potential 

Vqo = 7 ; (%i + x\ + xf) + bx i + cx 2 + dx 3 + e 
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is a non-degenerate second-order superintegrable system. The BD equations for the potential 
are 


(V n\ 


/I 0 0 0\ 


E22 


10 0 0 


(Vee\ 

E33 


10 0 0 


V 1 

Ul2 


0 0 0 0 


E 2 

E13 


0 0 0 0 


\vj 

\V,23/ 


\0 0 0 0/ 



Taking the origin as a generic point, the coefficients are related to k as a = V ee (0), b = V[i(0), 
c = P 2 (0), d = K 3 (0), and e = U(0). 

2.1 Equivalence classes 

The next step in the classification is to identify the appropriate equivalence classes. Clearly, we 
would like to consider two potentials as equivalent if they are related by a change of parameters. 
This is accomplished in the previous section by making a canonical choice of parameters as the 
value of u(x 0 ). We would also like systems to be equivalent if they are related by translations 
in the position variables; such a translation corresponds to moving the regular point. We would 
also like systems to be equivalent if they are related by, possibly complex, rotations. This con¬ 
dition will be studied extensively in this section. Finally, we would like the equivalence classes 
to include systems that are related via the Stackel transform or coupling constant metamor¬ 
phosis [10, 11, 17, 24], That is, suppose a Hamiltonian can be expressed as 77 = 77o + cdJ 
then the Stackel transformed Hamiltonians 77 = jjl-L will also be superintegrable, perhaps on 
a different conformally flat manifold. Under this transformation, classical trajectories and quan¬ 
tum wave functions as well as the corresponding symmetry algebras are essentially preserved, up 
to a change of parameters. Thus, we would like to consider two such systems as being equivalent. 
An analogous classification of equivalence classes for superintegrable systems in 2D has already 
been performed [20]. 

In order to consider Stackel equivalent systems, we focus our attention on conformally super¬ 
integrable systems [13]. For conformal superintegrable systems, the Hamiltonian is identically 0 
on trajectories. This can be accomplished mathematically by including the trivial added pa¬ 
rameter to the system. Thus, along a trajectory the Hamilton-Jacobi equation is expressed as 
77 — E = 0 but with 77 — 77 being the Hamiltonian. A function on phase space will then be 
a constant of the motion whenever its Poisson bracket is given by 

{A 77} = /(p,x)77. 

It is a direct computation to verify (see, e.g., [2, Lemma 4.1.4]) that a conformal integral of 
the motion A of 77 will be a conformal integral of the scaled Hamiltonian t/(x)77. Thus it is 
possible to transform any conformal superintegrable Hamiltonian on conformally flat space to 
a conformally superintegrable Hamiltonian on Euclidean space. Furthermore, since the action 
of the Stackel transform corresponds to multiplying the Hamiltonian by a function, it is clear 
that two Hamiltonians equivalent under this action will be equivalent conformal superintegrable 
systems, from the perspective of their symmetry algebras. 

Example 2. The Hamiltonian for the simple Harmonic oscillator p 2 — u; 2 x 2 corresponds to the 
following conformal Hamiltonian 


77 = p 2 - cu 2 x 2 - E. 
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Dividing by x 2 gives 


n 



E 


Clearly, trajectories that satisfy 7i = 0 will also satisfy 71 = 0, however the “energy” of the new 
Hamiltonian is given by the parameter co 2 . 


Conversely, any conformally superintegrable system that depends linearly on at least one arbi¬ 
trary constant can be transformed, see, e.g., [2, Theorem 4.1.8], into a superintegrable system by 
coupling constant metamorphosis. Therefore, classifying superintegrable systems on conformally 
flat manifolds up to Stackel transform is equivalent to classifying conformally superintegrable 
systems on Euclidean space. The determining equations for an integral of a conformal system 
are similar to (1) except that the additive constant is no longer trivial but instead essential for 
specifying the Hamiltonian. The corresponding structure equations are then 


/Vn\ 

H 2 2 

H 33 

V 12 

Hi 3 
\V,23/ 


A 


[Vee\ 

Hi 

H 2 

Hs 

V H / 


with 


(1 

-4S 1 - 

Rf - Rf 

2S 2 + R 12 

2 S 3 + R{ 3 

Af\ 

1 

2S 1 

+ R\ 2 

-AS 2 - R\ 2 - Rf 

2 S 3 + Rf 

Af 

1 

2S 1 

ARf 

2 S 2 + Rf 

-AS 3 - Rf - Rf 

Af 

0 

R{ 2 

-3 S 2 

Rf - 3S 1 

Q 123 

Af 

0 

R l 3 

-35 s 

Q 123 

Rf - 3S 1 

Af 

\0 

Q 123 

Rf - 3 S 3 

Rf - 3 S' 2 

Af) 


The functions Af are quadratic functions in the {Q, S, R} given in [2], 


2.2 Action of the conformal group 

So far, we have determined that superintegrable systems are uniquely determined by the value 
of the 10 functions {Q,R, S} at a regular point. We are interested in determining equivalence 
classes of systems that are stable under the action of the conformal group so we look at the 
induced action of this group on these 10 functions. Again, the exact formulas and derivations 
can be found in [2, 3] and here we state only the results necessary to understand the limits 
obtained in Section 3. 

A conformal change of variables can be generated by translations, which act trivially on the 
functions {Q, R, S} and inversions in the spheres of varying radii, which decompose the functions 
into a 7-dimensional representation {Q,R} and a 3-dimensional representation carried by {S}. 
Significantly, we can use a group transformation to set the values of {S} to any desired value, 
described in the following theorem. 

Theorem 1 (re-statement of [2, Theorem 4.2.18]). Given a conformally superintegrable system 
with function values {Qo, Ro, S 0 } at a regular point xo, there exists a conformal group motion 
that maps the function values to {Qo, Ro, 0} at the transformed regular point xq . 
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This theorem allows us restrict our attention to the action of the conformal group on the 
7-dimensional space {Q,R}. To understand the action, we consider the continuous generators 
of the conformal group, namely translations, scaling, rotations and Mobius transformations. 
As discussed before, translations act trivially on our representations. Scaling the coordinates 
corresponds to a scaling of the functions. More interesting is the action of rotations. To represent 
this action, we form weight vectors from the functions 

^±3 = R I 12 + \r 2 3 3 ± * (r? + \r1 3 ^ , Y ± 2 = - Rf) =F 2Q 123 ), 

Y±i = =f iRl 3 ), Y 0 = ~iVb(Rl 3 + Rf), 

so that the action of rotations are given by the standard raising and lowering operators for so(3). 
Using the isomorphism between 503 (C) and 3(2 (C), we obtain a covariant representation of the 
action via the polynomial 

9« = E (- 1)1 \[Q Y3 -Y^ ■ 

Here the action of SOs(C) is represented, via the isomorphism, as the standard action of SL 2 (C) 
as Mobius transformations 

^) = (ci + rf) e,(|±f). ( 2 ) 

Thus, the action of the conformal group on the functions {Q,R} can be represented by the 
action of GL 2 (C) = C x SL 2 (C) on the polynomials q(z ) via scaling and (2). Furthermore, 
the action of Mobius transformations on the coordinates can be represented locally as rotation 
combined with a scaling and so, if we find invariants that are closed under translation, scaling 
and rotations they will be automatically closed under the entire action of the conformal group. 
Finally, we have the following theorem. 

Theorem 2 (Theorem 4.2.21 of [2] in non-homogeneous coordinates). Given a conformal super- 
integrable system and a regular point, there is a local conformal transformation ( i.eexcluding 
translation of the regular point ) taking it to the regular point of another superintegrable system 
if and only if the roots of the corresponding covariant polynomial at the corresponding regular 
points are equivalent up to a general linear transform. 

Therefore, a key to understanding the classification of conformal superintegrable systems is 
to understand the classification of invariants of the roots of degree 6 polynomials. Furthermore, 
since we would like the equivalence classes to be stable under translation of the regular point, we 
require that any invariants obtained be closed under derivations. This significantly reduces the 
number of equivalence classes to 10 , exactly the expected number of equivalent superintegrable 
systems. For more details on the classification, we refer the reader to the original papers [2, 3]. 


3 Contractions 


For each of the contractions, we need only record the following data. The initial regular point 
xo, the rotation angles ti, t 2 , t%, the scaling parameter c and the final regular point yo- We 
note that the rotations can be represented in SL 2 (C) by the following matrices: 


p(Rs) 


' e it 3 /2 0 

0 e -*W 2 ’ 


p(R*) 


cos(f 2 /2) — sin(f 2 /2) 
sin(t 2 /2) cos(f 2 /2) ’ 
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p(Ri) 


cos(fi/ 2 ) — isin(ti/ 2 ) 

—*sin(ti/ 2 ) cos(ti/ 2 ) 


The scaling c can be encoded in the same manner by the matrix 


P(c) 


c " 1 / 6 0 

0 c -1/6 


In the following sections, unless otherwise mentioned, the angles are assumed to be set to 0 and 
the scale factor c to be set to 1. The action of the matrices on the polynomials are given by 


p(A) o q(z) = q(z), p{A) = 


as in ( 2 ), which transform the roots as 


da — b 

a = -. 

—ca + a 


To emphasize the dependence of the polynomials on the regular point, we write q(xi, x 2 , x^){z) 
where necessary. 

The action of the SL 2 (C) matrices can also be interpreted in terms of their action on a ste¬ 
reographic projection of the complex plane onto the unit sphere using the formula 


(X,Y,Z) 


2 x 


2 y 


1 + x 2 + y 2 ’ 1 + x 2 + y 2 ' 


-1 +x 2 + y 2 \ 
1 + x 2 + y 2 ) 


For real t\. the action of is to rotate the roots around the X axis by an angle t\ counter¬ 

clockwise (i.e., from Y to Z ). Similarly, the action of R 2 (t 2 ) is to rotate the roots around the Y 
axis clockwise (i.e., from X to Z). Finally, the action of 7 ? 3 (t 3 ) is to rotate the roots around 
the Z axis clockwise (i.e., from Y to X). 

The action of these rotations on M 3 can be recovered as 


Ri 


Rs 


10 0 
0 cos(fi) — sin(ti) , 
0 sin(fi) cos(ti) 

cos(fs) — sin ^ 3 ) 0 

sin(t3) cos(t3) 0 . 

0 0 1 


R 2 


cos(i2) 0 sin(f2) 
0 1 0 
— sin(f2) 0 cos(i2) 


In total, the change of coordinates for the potential is given by 


y = cR 3 {t 3 )R 2 (t 2 )Ri(ti)(-K - x 0 ) + y 0 . 


Notice that moving the regular point corresponds to translating the coordinates. 

For the remainder of this section, we give the limits between the systems. The first few case 
involve three or fewer roots, or (as in the case of [3111b]) four roots with a fixed cross ratio. The 
limits between these cases can be achieved solely with the action of GL 2 (C), that is, without 
appealing to translation the regular point. For most of these case the limits can be achieved by 
fixing one root and collapsing the rest together at infinity (or another prescribed point in C*). 
Several of the limits are elaborated to give a more complete understanding of the process. 



















J.J. Capel, J.M. Kress and S. Post 


3.1 [6] V A to [0] Vo 

The covariant polynomial for the [0] equivalence class is simply 0, so to contract down to this 
equivalence class, we need only scale by a factor of e _1 and move the regular point appropriately. 
For Va, the regular point is already xo = (0,0,0) so it is unaffected by the scaling. The 
contraction is then 


h=t 2 = t 3 = 0 , c = e\ x 0 = ( 0 , 0 , 0 ), y 0 = ( 0 , 0 , 0 ), 
q A ( 0 , 0 , 0 )(z) = -iz 6 -> q o ( 0 , 0 , 0 )(z) = 0 . 

We start with the potential 


V A = a Q (xi 2 + x 2 2 + x 3 2 ) + -^(aq - ix 2 ) 3 ^j + b (^c i + |(aq - ix 2 ) 2 ^j 


+ c[x 2 — -(xi — ix 2 ) 2 ) + dx 3 + e. 


The conformal change of coordinates y = e _1 x introduces a factor of e“ 2 to the second order 
terms in the Hamiltonian (due to the conformal scaling of the metric) and so we consider the 
rescaled potential V A = e 2 V A . Adjusting the parameters to be 

a = V A ,ee(yo) = € 4 V A ,ee(x o) = e 4 a, b = VA,i(yo) = e 3 Vk ) i(x 0 ) = e 3 b, 

C = V Aj 2 (y 0 ) = 6 3 Va, 2 (x 0 ) = e 3 c, d = VA, 3 (yo) = e 3 V) 4 , 3 (xo) = e 3 d, 

e = Vk(yo) = £ 2 Va(xo) = e 2 e, 
gives the potential 

V A = e 2 V A = d Q(yi 2 + y 2 2 + y 3 2 ) + ^(yi - +b(y x + |(yi - iy 2 ) 2 ) 

+ c fy 2 - j(yi - m) 2 ] + dy 3 + e, 


which tends in the e —» 0 limit to the potential 

Vo = | (y? + yl + yi) + &yi + %2 + %3 + e. 


3.2 [51] Vvn to [6] Va 


The regular points for Vyu and V A are the origins, xo = 0 and yo = 0. The change of variables 
is given by 


c = 9/2e 3 , t A = 0, t 2 = t 3 = i ln(e), 

which transforms the polynomial 

yvn(0, 0, 0)(z) = — 36iz -A q A (0, 0,0)(z) = — iz 6 . 

Consider first the roots of the (projective) sextic gyn(0,0,0) = —36 iz. There are five at 
z = oo and one at z = 0. 



Re(z) 


0 
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In terms of a stereographic projection of the complex plane z = x + iy (with the plane of 
projection intersecting the equator) the action of I? 2 ( 0 ) is to rotate the sphere around the y- 
axis (i.e., the imaginary axis which is directed into the plane) by the angle £ 2 - So consider 

P(^2(|))=ff f)- 

\V2 V2 / 

Under the action described above the covariant polynomial becomes 
p(R 2 (f )) 0 <?vii = “ !). 

which has brought the fifth root down from infinity down to z = — 1 and moved the root at zero 
up to z = 1. 


00 



Re(z) 


Setting tj, = i ln(e) gives the matrix 

/ e -1 / 2 0 

p(R 3 (iln(e))) = ^ Q el/2 
On the polynomial this induces the change 

(p(-R 2 (f))p(-R 3 (*lne))) oqvii = - ^ie~ 3 (z + e) 5 (z- e). 

So finally, scaling by c = §e -3 gives 

(cp(R 2 (^))p(R 3 {i In e))) o q YU = -i(z + e) 5 (z - e) -> q A = -iz % as e0. 

This limit can be visualized as all the points on the sphere (except 00 ) being drawn down towards 
the origin. 




Re(z) 


Or equivalently the roots can remain stationary while the sphere descends into the plane of 
projection. From this point of view, the image of any point on the sphere (expect 00 ) end 
up eventually being projected inside the circle where the plane and sphere intersect, which is 
shrinking to the origin in the limit. 
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3.3 [33] V 00 to [6] V A 

Here, we would like to take the two roots of goo, namely 0 and oo, and move them both to 0. 
This is going to be essentially the same contraction as the previous one. The contraction is 

c = 3/4e~ 3 , ti=iln(e), t 2 = *3 = 0. 

The regular points are 

x 0 = (0,0,1), y 0 = (0,0,0). 

The covariant polynomial is transformed as 

q O o(0,0, l)(z) = 6iz 3 -* q A (0,0,0)(z) =-iz 6 . 


3.4 [411] V IU to [51] Vvn 

The covariant polynomial for Vni is given by gni = —3(1 + 3 z 2 ). In order to transform this 
polynomial into gyn = —36*;? we need to take one of the finite roots and send it to zero and the 
other to infinity. The contraction that accomplishes this is 


^L e _ 2 
32\/3 ’ 


ti = --ir, t -2 = 7r, * 3 = —iln(e). 


The regular points for each system do not move and they are 


x 0 = (0,1,0), y 0 = (0,0,0). 


The limit of the polynomial is then 

<jfm(0,1,0)(z) = — 9z 2 — 3 —>• gvn(0,0,0)(z) = — 36*2? 


3.5 [3111b] Vvi [51] Vvn 

The necessary contraction is given by 

c = Yg e_1 > tl = \' *2 = 0, f 3 = —^ln(e), 

x 0 = (0,0,2*), y 0 = (0,0,0). 

In this limit, all of the roots coalesce to infinity except for a single root at zero. The limit of 
the polynomial is then 

</vi(0,0, 2*)(z) = 3iz 6 + 3z 3 -> 9vn(0, 0, 0)(z) = — 36iz. 


3.6 [3111b] Vvi to [33] V 00 

For this contraction, we would like to move the roots of the polynomial qyi by sending the three 
non-zero roots to infinity and leaving the 0 root fixed. This can be accomplished by choosing 
an imaginary, singular value for 1 3 which has the effect of scaling the variable z. The required 
contraction is then 

* 3 = — *ln(e), c= \’ x 0 = (0,0,2), y 0 = (0,0,1). 

The polynomials transform as 

<7vi(0,0, 2)(z) = 3*2: 3 (l + z 3 ) -* g O o(0,0, l)(z) = 6 iz 3 . 
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3.7 [3111a] V u to [411] F v 

From this case onward the polynomials involve four or more roots and require additional infor¬ 
mation about the cross ratios of the roots in order to differentiate between equivalence classes. 
Translation of the regular point allows the cross ratios of these roots to change. Without transla¬ 
tion of the regular point the only way singular limits could be taken is by collapsing all, save one, 
of the roots together. This means that, starting at a point which has represents root multiplic¬ 
ities [3111] the only possible limits without translation of the regular point are [51],[33] or their 
degenerations. In all of the examples that follow, translation of the regular point will be neces¬ 
sary to modify the cross ratios (i.e., move the roots in a way which the Mobius transformations 
cannot) and, in effect, rescue some of the roots in the singular limit. 

To contract the covariant polynomial of Vn, we observe that the polynomial evaluated at 
(xi,X 2 ,x 3 ) is 


( 2 3 

—H--v—+ 

X 3 X\ + 1X2 

Specializing the first two coordinates to (x\,X 2 ) = 


X] - ix 2 \ 

(aq + ix 2 ) 2 J 

(0, — i ) gives 


<?n(0, -i,x 3 )(z) = 3 i 


X3 


z J + 3z 2 - 1 


and so by moving the regular point to 00 , we obtain the desired contraction. The regular points 
and scaling parameter is then chosen as 


c = i, x 0 = (0, —i, e 1 ), y 0 = (0,1,0). 


This contraction gives the required limits of the polynomial, 

q u {0,-i,l){z) = 3i(2z 3 + 3z 2 -1) -> <?v(0,1, 0)(z) = 9z 2 - 3. 


3.8 [3111a] V u to [3111b] V V i 

There are two distinct equivalence classes of systems whose covariant polynomials have the root 
structure [311]; these equivalence classes are distinguished by the characteristic that the cross 
ratio of the root be exp(i7r/3) (case [3111b]) or not (case [311a]) [2, 3]. For the limit between 
these two systems, we begin with a polynomial with a triple root and three simple roots with 
a cross ratio of the roots away from 0, 1, 00 or exp(±i7r/3). We would like to find a limit 
that will contract this polynomial to one with the same root structure but with a cross ratio of 
exp(*7r/3). We begin with the covariant polynomial for V\\ which has a triple root at infinity 
and move this root to 0 by rotating around the X 2 axis by it. We then move the x\ coordinate 
of the regular point to infinity and perform a singular rotation about x 3 via 

L 2 = 7L t 3 = ^hi(e), x 0 = (e _1 ,0,3e-2), y 0 = (0,0,2). 

This limit transforms the polynomials as (from e = 1 to e = 0) 

gii(l,0,l)(z) = 3i(l + 3z 2 + 2z 3 ) -> qvi(0,0,2)(z) = 3iz 3 (z 3 + l). 

3.9 [111111c] V lv to [411] Vy 

For this contraction, we take the equivalence class described by a polynomial with 6 simple roots 
satisfying the multi-ratio condition, a discriminating condition for the [111111] root structu¬ 
re [2, 3], to one with two simple finite roots and the remaining at infinity. The coordinate limit 
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is given by 

— 1 7T 

c=—, ii = 0, t 2 = -~, t 3 = -i ln(e). 

x 0 = > yo = (0,1,0). 

The covariant polynomial contracts as 

qiv(xi,x 2 ,x 3 )(z) = + ^ ->• gv(0,1,0)(z) = 9z 2 — 3. 

x 3 4x 2 


3.10 [111111c] Viv to [3111b] Vvi 


As seen above, the covariant polynomial for the potential Vjv is given by 

^ Q 

qiv(xi,x 2 ,x 3 )(z) = —z 3 + -— (z 2 + 1)' • 

X 3 4X2 

An obvious easy way to achieve a triple root is to let the second coordinate x 2 tend to infinity. 
However this has the unintended side effect of collapsing the other three roots at infinity. The 
singular rotation R 3 (i\n(e)) can be used to counteract this growth, provided the parameter e 
correctly balances the growth of x 2 . This action only serves to accelerate the roots tending to 
zero, and hence gives us the limit we require. 

This contraction is accomplished by a singular rotation in around the X3-axis and a singular 
translation of the regular point. The contraction is 


t 3 = i In e, 


x ° t 4 e 3 ’ 2 


yo = (0,0,2). 


The covariant polynomials contract as 


s(V 2 _i_ 11 3 

qw{xi,x 2 ,x 3 ){z) = —-1-—- ->• qvi (0, 0, 2) (z) = 3iz 3 (z 3 + l). 


2:3 


4x 2 


3.11 [111111b] Vi to [3111a] V u 


Here we would like to send 3 of the finite roots to infinity. The covariant polynomial for V\ is 


qi(x 1 ,x 2 ,x 3 )(z) = 


3(xi+ix 2 )z 6 9 (xi — ix 2 ) z 4 6 iz 3 9 (xi + ix 2 ) z 2 3(xi — ix 2 ) 


4xix 2 


+ 


4xiX2 


+ 


3?' / 2 i \ 3 3 / 2 | i \ 3 | 6* 

4^ (z - 1 ) + 4a?2 2 +1 > 


®3 

3 


+ 


4xiX2 


+ 


4xix 2 


Rotations around the z axis have the effect of multiplying the variable z by e~ 2its and the entire 
polynomial by a factor of e 6 ** 3 . Thus if we take a purely imaginary angle, the z 3 variable will 
be unchanged and it will be possible to send the coefficient of z 4 and z e to 0. In order to retain 
the lower order terms, it is necessary to move the regular point so as to offset the singular limit. 
The required contraction is 


^,ij, yo = (i,o,i). 

The polynomials transform as 

<7i(xi,x 2 ,x 3 )(z) = - l) 3 + ^( z2 + !) 3 + 

q n (l,0,l)(z) =3i(l + 2z s + 3z 2 ). 


t 3 = —i lne, 


x 0 = 


e 2 + 1 iie 2 - 


2e 


2e 
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3.11.1 “Geometric inspired limit” 

This contraction can be carried out in a manner similar to the previous one. By moving the 
regular point correctly we can collapse three roots at z = 0 and three roots are z = oo. Using 
the same singular rotation as before, the collapse of one of these triplets of roots can halted, 
and the other accelerated. 

The covariant polynomial for V\ is given by 

qi(x 1 ,x 2 ,x 3 )(z) = ^-(z 2 - l ) 3 + ^(- 2 + 1 ) 3 + ^ 3 - 

By choosing to work on the hypersurface X 2 = —ix\ we ensure that q\ has a double root at 
z = 0. The choice x\ = 77 , x = —ir] gives 

• w , 3i 2 ( zA 4 z 3\ 

<n{v, -iv, x 3 ) (z) = —z —1 - 1 — ■ 

2 \rj x 3 r)J 

Sending rj to infinity gives the limit collapses a single root to join the double root at zero, and 
sends the other three roots to infinity. Using a singular rotation around the x 3 axis the root 
tending to zero can be saved, and the roots tending infinity accelerated. 

The required contraction is therefore 

c= l> t 3 = -ilne, x 0 = (e _1 ,-ze _1 ,4), y 0 = Q,^,2). 

The polynomials transform as 

qi(xi,x 2 ,x 3 )(z) = ^(, 2 - 1 ) 3 + i ^(, 2 + 1 ) 3 + |, 3 ->• 

, 2 ^ (z) = 3 iz 3 + 9 iz 2 . 

3.12 [111111b] V l to [111111c] Viv 

To obtain this contraction, we move the regular point and scale the potential. The required 
contraction is 

c=~, x 0 = (1, e, e), yo = (0,1,1). 
e 

This transforms the covariant polynomials as 

qi(xi,x 2 ,x 3 )(z) = ^(, 2 - 1 ) 3 + ^(, 2 + 1 ) 3 + |, 3 ->• 

fflv( 0 , 1,1 )(z) = 6iz 3 + ^(z 2 + l) 3 . 

3.13 [111111a] V S w to [111111b] Vi 

The case [111111a] covers the most generic forms for the covariant polynomials. A particularly 
nice example of the polynomials in this class is 

9 

<?sw(x o )(20 = -z(z - l){z + 1 ){z - i)(z + i), x 0 = ( 1 , 1 , 1 ). 

However, this polynomial doesn’t show the full regular point dependence of qsw(z) 
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One possible contraction is given by a scaling by e 1 and translation of the regular point to 
the origin 

c=~, x 0 = (e, e, e), y 0 = ( 1 , 1 , 1 ). 


Following the regular point (e, e, e) the covariant polynomial for the is given by 

3?'(1 — e 4 )(l + i) 36 e A z 9(1 — e 4 )(l + i)z 2 6 i(l — e 4 )z 3 

q S w(e,e,e)(z) = ^ + (] - 9 e 4 ) 4(1 - 9 e 4 ) (1 - 9 e 4 ) 

9*(1 — e 4 )(l + i)z A 36 e 4 z 5 3(1 - e 4 )(l + i)z 6 


+ 


4(1 - 9e 4 ) 


(1 - 9e 4 ) 4(1 - 9e 4 ) 


Here e must follow a path from 1 to 0 in the complex plane, avoiding points where e 4 = -. 
The limit of this contraction transforms the covariant polynomial to 

3i / 9 >3 3/o \ 3 


q I (l,l,l)(z) = —( z2 - 1) +-{z 2 + l) + 6 iz 3 . 


4 Conclusion 

By using the invariant classification of conformally superintegrable systems on conformally flat 
3D manifolds, we have been able to give explicitly the coordinate limits between the 10 equiv¬ 
alence classes. We have therefore proven that all such systems are limits of the generic sys¬ 
tem, Vsw- We would like to point out that this analysis relied on the potentials begin complex, 
beginning with the fact that we are determining only equivalence classes under the Stackel 
transform which does not necessarily preserve real forms and continuing into the classification 
using the action of the conformal group. It would therefore be interesting to understand the 
classification of the real forms of these potentials and their corresponding limits, which would 
be non-trivial applications of the complex theory. 

One of the most interesting avenues of this research is the connection between superintegrable 
systems, orthogonal polynomials (OPs) and special functions. From the beginning of the study 
of superintegrable systems, going back to the harmonic oscillator, Kepler Coulomb system and 
the discovery in ’65 of multiseparable systems by Smorodinsky, Winternitz and collaborators [4, 
21, 27], there has been a clear connection between superintegrable systems and the special 
functions of mathematical physics. In particular, the wave functions of superintegrable systems 
have been expressed in terms of OPs [15] and all superintegrable systems are conjectured to have 
this exactly-solvable nature [26]. More recently, this connection has been extended by studying 
the representation theory of the symmetry algebras associated with superintegrable systems and 
interbasis expansion coefficients [5, 6 , 7, 8 , 9, 16, 18, 25]. In [19], the Askey scheme of classical 
hypergeometric OPs in one variable as well as their limits were related to superintegrable systems 
in 2D and the contractions of their symmetry algebras. Work is currently ongoing to apply the 
limits obtained in this paper to find contractions of the associated quadratic algebras and their 
polynomial representations. 

Finally, we would like to mention that as is evident from the classification theory, the su¬ 
perintegrable systems are uniquely determined by the functionally linearly independent set of 
second-order Killing tensors. Therefore we expect, as in the 2D case [14], that the quadratic 
algebras and their contractions will be completely determined by contractions of the underlying 
Lie algebra. However, this is still an open question and one which we anticipate will be aided by 
understanding the connection between the invariant classification and the quadratic algebras. 
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Figure 1. Limiting diagram. 


A Potentials 


Here we list the potentials discussed above. Each of the limits above has been chosen so that the 
limit potential is in the form given below. Note that the parameters are chosen is for aesthetic 
reasons, to give the potentials in the most simple form. They are related by a linear transforma¬ 
tion to the parameters {H(xo), bi(xo), b^xo), V^xo), F ee (xo)}, which of course depend on the 
regular point xo- 

Note that by is Stackel equivalent to bin by choice of the Stackel multiplier U = (xi + ix 2 ) -2 : 


b 0=2 ( x ’i + x 2 + x l) + + cx 2 + dx 3 + e, 

Va = a Q (xi 2 + x 2 2 + xs 2 ) + (xi - ix 2 ) 3 ^ + b ^zq + |(xi - ix 2 ) 2 ^ 
+ c ^x 2 — |(xi — ix 2 ) 2 ^ + dx 3 + e, 

Vqo = o(4x 2 + 4x| + x 2 ) + bx\ + cx 2 -\—^ T e, 
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Vyu = a(x i + ix 2 ) + b (3(xi + ix 2 ) 2 + x 3 ) 

+ c (l6(xi + ix 2 ) 3 + (x 3 - ix 2 ) + 12x 3 (xi + ix 2 )) 

+ d(5(xi + ix 2 ) 4 + x 2 + xl + x§) + 6((xi + ix 2 ) 2 x 3 ) + e, 

Vyi = a(x\ — 2(xi — ix 2 ) 3 + 4(x 2 + x^)) + b(2xi + 2 ix 2 — 3(x 3 — ix 2 ) 2 ) 


+ c(xi - zx 2 ) + + e, 


Vy — o(xj + x 2 + 4x 2 ) + &x 3 + 


, Xl - ix 2 

+ d -:-ry + e, 


(. xi+ix 2 ) 2 (xi + zx 2 ) 3 


Viv — a(4x 2 + x 2 + X 3 ) + fcxi 4 —2 4 —9 4~ 6 , 


XS Xq 


Km = a(x 2 + x| 4- x§) 4- 


cx 3 


1 /-v*2 _ q™2 

,0/1 "T 0/0 OJbo 

+ d-h - 2 ,/ +e, 


(xi 4- ix 2 ) 2 (xi 4- ix 2 ) 3 (xi 4- ix 2 ) 4 


Vn = a(x 1 2 + x 2 2 4- x 3 2 ) + b ^ Xl iX2 } Q + c 


1 


, 1 

4" d—~ 4“ e, 


Ki = a(x 2 4- X 2 4- x|) 4--^4-—^4--^4-e 


(xi 4- tx 2 ) 3 (xi 4- ix 2 ) 2 x 3 2 

c 

ry> ^ nf* 

Jb OO) O/ 


Ksw = 


c 


ri 


2 + 2 + 2 4- 2 + 


(l 4- xf 4- xf 4- xf) x \ x 2 x 3 (—1 4- xf 4- xf + xf) 


2 • 
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